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Abstract. Massive black holes are nowadays believed to reside in most local galaxies. Studies 
have also established a number of relations between the MBH mass and properties of the host 
galaxy such as bulge mass and velocity dispersion. These results suggest that central MBHs, 
while much less massive than the host (~ 0.1%), are linked to the evolution of galactic structure. 
When did it all start? In hierarchical cosmologies, a single big galaxy today can be traced back 
to the stage when it was split up in hundreds of smaller components. Did MBH seeds form 
with the same efficiency in small proto-galaxies, or did their formation had to await the buildup 
of substantial galaxies with deeper potential wells? I briefly review here some of the physical 
processes that are conducive to the evolution of the massive black hole population. I will discuss 
black hole formation processes for 'seed' black holes that are likely to place at early cosmic 
epochs, and possible observational tests of these scenarios. 

Keywords, cosmology: theory, galaxies: formation, gravitation, black hole physics, quasars: 
general 



1. Introduction 

I will focus here on the formation and evolution of massive black holes (MBHs) , in high- 
rcdshift galaxies, and their symbiotic evolution with their hosts. MBHs weighing million 
solar masses and above have been recognized as the engines that power quasars detected 
at early cosmic times. Dynamical evidence also indicates that MBHs with masses in the 
range Mbh ~ 10 6 — 1 9 Mq ordinarily dwell in the centers of most nearby galaxies 
( Ferrarese fc Fordl 2005). MBHs populate galaxy centers today, and shone as quasars 
in the past; the quiescent MBHs that we detect now in nearby bulges would be the 
dormant remnants of this fiery past, suggesting a single mechanism for assembling MBHs 
and forming galaxies. The surpri singly clear correlatio ns between MBH masses and the 
properties of their host galaxies (|Gultekin et al. 2009, and references therein) may well 
extend down to the smallest masses, for example, the dwarf Seyfert 1 gala xies POX 52 
and NGC 4395 are thoug ht to contain a MBH of mass M B h ~ 10 5 M Q (jBarth et al. 
120041: iPeterson et al.ll2005r ). At the other end, however, some powerful qu asars, powered 
by super-massive black holes with masses ~ 10 9 M Q jBarth et all 120031. have already 
been detected at z > 6, corresponding to a time less than a tenth of the age of the 
Universe, roughly one billion years after the Big Bang. 



We are therefore left with the task of explaining the presence of billion solar masses 
MBHs when the Universe is less than 1 Gyr old, and of much smaller MBHs lurking in 
13 Gyr old galaxies. The outstanding questions concern how and when "seed" MBHs 
formed, the frequency of MBHs in galaxies, and how efficiently MBH seeds grew in mass 
during the first few billion years of the i r lives . 

The "flow chart" presented by Reesl ( 1978 ) still stands as a guideline for the possible 
paths leading to formation of massive MBH seeds in the center of galactic structures. In 
this paper I will concentrate on the currently favoured scenarios. 
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2. Massive black hole 'seeds' 



2.1. Population III remnants 

One of the most popular scenarios for MBH formation associates MBH seeds with the 
remnants of the first generation of stars, formed out of zero metallicity gas. The first stars 
are expected to form in minihalos, M m ; n w 10 6 M© collapsing at z ~ 20 — 50 from the 
highest peaks o f the primordial de nsity field, where cooling is possible by means of molec- 
ular hy drogen (iPalla et alJl2002l ) . Simulations of the collapse of prim ordial molecular 



clouds (jBromm et al.lll999L Eool lAbel et ai]l200d : FYbshida et al.ll2006l ) suggest that the 



first generation of stars contained many 'very massive stars' (VMSs) with > 100 M Q . 
This is because of the slow subsonic contraction set up by molecular hydrogen cooling. 
Further fragmentation into sub-components is not see n although it is not clear if this is 
a numerical effect, rather than due to the gas physics ( Glover et alJl2008h . 

If the first stars retain their high mass until death, they will collapse after a short (w 
Myrs) life-time. The final fate depends on the exact mass of the star. Between 25 and 140 
M , low- metallicity stars form black holes directly. It is likely, however, that such light 
black holes wander within their hosts, dynamically interacting with stars of similar mass, 
rather than settle at the center of the galaxy's potential well. Between approximately 140 
and 260 M lies the domain of pair instability supernovae. The se objects are completel y 
disrupted by nuclear-powered explosions, leaving no remnants ( Kudritzki fc Puis] 2000h . 
Going to still more massive stars (over 260 Mq) the nuclear energy released when the star 
collapses on the pair instability is not sufficient to reverse the implosion. A m assive black 



hole, containing at leas t half of the initial stellar mass, is born inside the star (jBond et al 



1984; Fryer et al.ll200lh . A numerous population of MBHs may have been the endproduct 
of the first episode of pregalactic star formation; since they form in density high peaks, 
relic MBHs with mass > 150 M would be predicted to cluster in the cores of more 
massive halos formed by subsequent mergers (jMadau fc Reedl200lh . 

Although this path to MBH formation seems very natural, large uncertainties exist on 
the final mass of PopIII stars. If the first stars form ed in multiples per halo the mass 
function would be less top-heavy (jClark et al.ll2008af ). And, if single objects are formed, 
the characteristic final mass at the end of the accretion process ca n be much less than 
the initial clump mass due to feedback effects ( McKee fe Tan 20081 ). 



2.2. Gas- dynamical processes 

Another family of models for MBH forma tion relies on the collapse of supermassive 
objects formed directly out of dense gas (e.g., Bromm fe LoebfeoOStlBegelman et al . 2006; 



Lodato fe Nataraian 20061 ). The physical conditions (density, gas content) in the inner 



regions of mainly gaseous proto-galaxies make these locii natural candidates, because the 
very first proto-galaxies were by definition metal-free, or at the very least very metal- 
poor. Enriched halos have a more efficient cooling, which in turn favors fragmentation 
and star formation over the efficient collection of gas conducive to MBH formation. It 
has been suggested that efficient gas collapse probably occurs only in massive halos with 
virial temper atures T v j r > 10 4 K un der metal-free conditions where the formation of H2 
is inhibited (jBromm fc Loeb|l2003l). or for gas en riched below the critical metallicity 



threshold for fragmentation (jSantoro fc Shulll 120061 ) . For these systems the tenuous gas 
cools down by atomic hydrogen only until it reaches T gas ~ 4000 K. At this point the 
cooling function of atomic hydrogen drops by a few orders of magnitude, and contraction 
proceeds nearly adiabatically. Finally, it has been recently suggested that highly turbulent 
systems are also likely to experience a limited amount of fragmentation, suggesting that 
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efficient gas collapse could pro ceed also in metal-enriched galaxies at later cosmic epochs 
( Begehnan fc Shlosmanll2009l ). 

In such halos where fragmentation is suppressed, and cooling proceeds gradually, the 
gaseous component can cool and contract until rotational support halts the collapse. 
In the most common situations, rotational support can halt the collapse before densi- 
ties required for MBH formation are reached. Additional mechanisms inducing transport 
of angular momentum are needed to further condense the gas until conditions foster- 
ing MBH formation are achieved. An appealing route to efficient angular momentum 
shedding is by global dynamical instabilities, such as the "bars-within-bars" mechanism 
( Shlosman et al.lll989t feegelman et al.ll2006l ). A bar can transport angular momentum 
outward on a dynamical timescale via gravitational and hydrodynamical torques, allow- 
ing the radius to shrink. Provided that the gas is able to cool, this shrinkage leads to 
even greater instability, on shorter timescales, and the process cascades. This mechanism 
is a very attractive candidate for collecting gas in the centers of halos, because it works 
on a dynamical time and can operate over many decades of radius. 

An alternative way of describing mass-inflow is via the Toomre stability parameter, 
Q. When Q approaches unity a disc is subject to gravitational instabilities. If the desta- 
bilization of the system is not too viole nt, instabilities lead to ma ss infall instead of 
fragmentation and global star formation ( Lodato fc Nataraian 20061 ). This is the case if 



the inflow rate is below a critical threshold M„ 



2a c 



the disk is able to sustain 



J-max — *<<-*c q 

(a c ~ 0.12 describes the viscosity) and molecular and metal cooling are not important. 
Such an unstable disc develops non-axisymmetric spiral structures, which effectively re- 
distribute angular momentum, causing mass inflow. This process stops when the amount 
of mass transported to the center is enough to make the disc marginally stable. 

After gas has efficiently accumulated in the center, what happens next? The gas made 
available in the central compact region can then form a central massive object. Depending 
on how fast and efficiently the mass accumulation proceeds, the exact outcome would 
differ. It is ge nerally thought that VMSs supported by radiation pressure evolve as a n n 
3 polytrope ( Hovle fc Fowler! [l963t Baumgarte fc Shapiro! 19991 : 



Saiio et al.ll2002h . The 



fate of a marginally unstabl e, maximally rotating VM S has been investigated numerically 
in full general relativity bv lShibata fc Shapiro! ( 2002 ). They found that the final object 
is a Kerr-like black hole (spin parameter «0.75) containing 90% of the stellar mass. The 
fate of an isolated VMS is therefore the formation of a MBH. 

If the mass accumulation is fast, however, the outer layers of VMSs are not the rmally 
relaxed during much of the main sequence lifetime of the star (jBegelman 2009). Such 
VMSs can have a complex structure with a convective (polytropic) core surrounded 
by a convectively stable envelope. After exhausting its hydrogen, the core of a VMS will 
contract and heat up until it suffers catastrophic neutrino losses and collapses. The initial 
black hole, with mass of a few M , formed as a result of core-collapse subsequently grows 
via accretion from the bloated envelopes that result ('quasistars', an initially low-mass 
black hole embedded within a massive, radiation-pressure-supported envelope; see also 



Begelman et al.l 12006! ) . Over time, the black hole grows at the expense of the envelope, 



until finally the growing luminosity succeeds in unbinding the envelope and the seed MBH 
is unveiled. The key feature of this scenario is that while the black hole is embedded within 
the envelope, its growth is limited by the Eddington limit for the whole quasistar, rather 
than that appropriate for the black hole mass itself. Very rapid growth can then occur 
at early times, when the envelope mass greatly exceeds the black hole mass. 

The masses of the seeds predicted by different models of gas infall and VMS structure 
vary, but they are typically in the range Mbh ~ 10 4 — 10 5 M (Figure [1]). 
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2.3. Stellar- dynamical processes 

Efficient gas collapse, leading to MBH seed formation, is mutually exclusive with star 
formation, as competition for the gas supply limits the mass available. Formation of 
'normal' low-mass stars, however, opens up a new scenario for MBH seed formation, if 
stellar-dynamical rather than gas-dynamical processes are at play. This first episode of 
efficient star formation can foster the formation of very com pact nuclear star clusters 
( Schneider et al.ll2006t Iciark et allfeoOSli lOmukai et alj[200^) where star collisions can 



lead to the formati on of a VMS, possibly leaving a MBH remnant with mass in the range 
~ 10 2 — 10 4 M©. iGaburov et al investigate the process of runaway collisions 

directly with hydrodynamical simulations and show that during collisions large mass 
losses are likely, at metallicities larger than 10 -3 solar. 

The growth of a VMS should be much more efficient at low metallicity. Low metallicity 
can modify the picture in different ways. First, at sub-solar (but non-zero) metallicity, 
stars with masses > 40M© are thought to collapse directly into a black hole without 
exploding as supernovae. Second, the mass loss due to winds is much more reduced in 
metal-poor stars, which greatly helps in increasing the mass of the final remnant. 

Devecchi fc Volonteri ( 20091 ) investigate the formation of MBHs, remnants of VMS 



formed via stellar collisions in the very first stellar clusters at early cosmic times. The 
main features of the model are as follows: we consider halos with virial temperatures 
T v i r >10 K after the first episode of star formation, hence with a low, but non-zero, 
metallicity gas content. This set of assumptions ensures that (i) atomic hydrogen cooling 
can contribute to the gas cooling process, (ii) a UV field has been created by the first 
stars, and (iii) the gas inside the halo has been mildly polluted by the first metals. The 
second condition implies that at low density H2 is dissociated and does not contribute to 
cooling. The third condition ensures that gas can fragment and form low-mass stars only 
if t he gas density is above a certain threshold, n crit: z (which depends on the metallicity, 



Santoro fe Shull 20061 ): this causes only the highest density regions of a proto-galactic 



sec 

disc to be prone to star formation. 

In Toomre- unstable proto-galactic discs, such as those described in the previous section, 
instabilities lead to mass infall instead of fragmentation into bound clumps and global 
star formation in the entire disk. The gas inflow increases the central density, and within 
a certain, compact, region n > n cr i t .z- Here star formation ensues and a dense star 
cluster is formed. At metallicities ~ 10~ 4 — 10~ 5 solar, the typical star cluster masses are 
of order 10 5 M© and the typical half mass radii ~ 1 pc. Most star cluster ther efore go 



into c ore collapse in < 3 Myr, and runaway collisions of stars form a VMS (see IShapiro 



20041 ) . leading to a MBH remnant. 

As the metallicity of the Universe increases, the critical density for fragmentation 
decreases and stars start to form in the entire protogalactic disk. As gas is consumed by 
star formation, the inflow of gas is no longer efficient, more extended clusters form, and 
the c ore collapse timescale increas es until it exceeds the main sequence lifetime of massive 



stars . iDevecchi fc Volonteril (2009) find that typically a fraction ~ 0.05 of protogalaxies at 



z ~ 10—20 form black hole seeds, with masses ~ 1000—2000 M©, lea ding to a mass density 



in see ds of a few ~ 10 2 M©/Mpc 3 . Most of the assumptions in [Devecchi fc Volonteri 



(2009) have been conservative, but still the population of se eds is comparabl e to th e case 



of Population III star remnants discussed, for instance, in IVolonteri et al.l (|2003l ). The 
fraction of high-redshift galaxies seeded with a MBH is about a factor of 10 below the 
direct collapse case presented in Volonteri et al.l ( 2008 ). where a seed was assumed to 



form with a 100% efficiency whenever a protogalaxy disk was Toomre unstable. 
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3. Observational signatures of massive black hole seeds 

What are the possible observational tests of MBH formation scenarios? Figured] shows 



three mass functions for three different MBH 'seed' mod els: direct collapse (IVolonteri et al 
2008 ) , runaway stellar me rgers in high-redshift clusters ( Devecchi fc Volonterill2009h . and 



Population III remnants (jVolonteri et al.ll2003l ). These are the initial conditions that we 
would like to probe. In this section, I will focus on discriminants for two different, ex- 
treme, scenarios: 'light seeds', derived from Population III remnants, and 'heavy seeds', 
derived from direct gas collapse. 




lg(M BH /M a ) 

Figure 1. Mass function of see d MBHs for three different formation scenarios: di- 
rect collapse (|Volonteri et all 120081 . left), runaway stellar merg ers in high-redshift clusters 
(Devecchi & Volontcri 2009, center), and Population III remnants (|Madau fe R ccs 2001, right). 
Note the different y-axis scale for the Population III case. 



3.1. Tracing MBHs at the earliest times 

Sinc e during the quasar e poch (z w 3 — 4) MBHs increase their mass by a large fac- 
tor (|Yu fc Tremainell2002r ). signatures of the seed formation mechanisms are likely more 
evident at earlie r epochs. If accretion is e fficient, it washes out any trace of the initial 
seed population. Volonteri fc Gnedin ( 20091 ) find that of the contribution to the quasar 
luminosity budget at z > 6 is dominated by MBHs with mass < 10 6 M Q . Such small, 
low-luminosity MBHs do not contribute to the bright end of the luminosity function of 
quasars, and are therefore difficult to account from simple extrapolations of the luminos- 
ity function of quasars. These small holes are not hosted in extremely massive galaxies 
residing in the highest density peaks (5 to 6-cr peaks), but are instead found in more 
common, "normal" systems, ~ 3 — a, peaks. Future generation of space-based telescopes, 
such as J WST and IXO, are likely to detect and constrain the evolution of the population 
of accreting massive black holes at early times (z>10). 
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3.2. Gravitational 



Detec tion of gravitational waves from seeds merging at the redshift of formation (jSesana et al 



20071 ) is probably one of the best ways to discriminate among formation mechanisms 



LISA in principle is sensible to gravitational waves from binary MBHs with masses in 
the range 10 3 — 10 6 M Q basically at any redshift of interest. A large fraction of coa- 
lescences will be directly observable by LISA, and on the basis of the detection rate, 
c onstraints can be pu t on t he MBH formation process. 

Sesana et al. (|2007l ) and lArun et all (|2009h analyze merger hist o ries of MBHs along 



besana et al.l (|^UU i\ ) and lArun et al.l l|zt)t)!j| ) analyze merger hist o ries ot Metis along 
the hierarchical build-up of cosmic structures ( Volonteri et al.|[2003h . Sesana et alj ( 2007t ) 



find that a decisive diagnostic is provided by the distribution of the mass ratios in binary 
coalescences. Models where seeds start large predict that most of the detected events 
involve equal mass binaries. A fraction of observable coalescences, in fact, involve MBHs 
at z > 10, when MBHs had no time to accrete much mass yet. As most seeds form 
with similar mass, mergers at early times involve MBH binaries with mass ratio ~ 1. In 
scenarios based on Population III remnants, z > 10 mergers involve MBHs with mass 
below the LISA threshold. The detectable events happen at later times, when MBHs have 
already experienced a great deal of mass growth yielding a mass ratio distribu tion which 
is flat or features a broad peak at mass ratios ~ 0.1 — 0. 2. lArun et all (|2009h determine 



the detectability of events by using a code that includes both spin precession and higher 
harmonics in the gravitational-wave signal, and carrying out Monte Carlo simulations 
to determine the number of events that can be detected and accurately localized in 
these population models. LISA will detect a significant fraction of all mergers: almost 
all mergers will be detected in heavy-seed scenarios, and nearly half of all mergers in 
small-seed scenarios. For heavy-seed models most mergers detectable with high signal- 
noise-ratio (S/N>10) occur in the redshift range 3<z<8, with a peak around z ~ 5. In 
the case of light seeds, mergers are roughly uniform in z over the range 4<z<10. 

3.3. Massive black holes in low-mass galaxies: occupation fractions and Mbh — o" 

The repercussions of different initial efficiencies for seed formation for the overall evo- 
lution of the MBH population stretch from high-redshift to the local Universe. The 
formation of seeds in a ACDM scenario follows the cosmological bias. As a consequence, 
the progenitors of massi ve galaxies (or cluster s of galaxies) have a higher probability of 
hosting MBH seeds (cfr. iMadau fc Reesl 2001 ). In the case of low-bias systems, such as 



isolated low-mass galaxies, very few of the high-z progenitors have the deep potential 
wells needed for gas retention and cooling, a prerequisite for MBH formation. 

The signature of the efficiency of the formation of MBH seeds will consequently be 
stronger in low-mass galaxies. Figure [2] (bottom panel) shows a comparison between the 
observed Mbh — o relation and the one predicted by different models (shown with c i rcles) , 



and in particular, from left to right, two models based on the lLodato fc Nataraianl (|2006l ) 
seed masses, and a third model based on lower-mass Population III star seeds. The upper 
panel of Figure [2] shows the fraction of galaxies that do host a massive black hole for 
different velocity dispersion bins. This shows that the fraction of galaxies with a MBH 
increases with increasing halo masses at z — 0. A larger fraction of low mass halos are 
devoid of central black holes for lower seed formation efficiencies. This is one of the key 
discriminants between different 'seed' scenarios. 

Simple arguments lead us to believe that MBHs might inhabit also the nuclei of dwarf 
galaxies, such as the satellites of the Milky Way and Andromeda, today (Van Wassenhove 
et al. 2009). Indeed, one of the best diagnostics of 'seed' formation mechanisms would 
be to measure the masses of MBHs in dwarf galaxies. As MBHs grow from lower-mass 
seeds, it is natural to expect that a leftover population of progenitor MBHs should also 
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Figure 2. The Mbh— velocity dispersion (<r c ) relation at z = 0. Every circle represents the 
central MBH in a halo of given er c . Observational data are marked by their quoted errorbars, 
both in (7 C , and in A/bh (Tremaine et al. 2002). Left to right panels: Q c — 1.5, Q c ~ 2, Population 
III star seeds. Top panels: fraction of galaxies at a given velocity dispersion which do host a 
central MBH. 



exist in the present universe. As discussed above, the progenitors of massive galaxies 
have a high probability that the central MBH is not "pristine" , that is, it has increased 
its mass by accretion, or it has experienced mergers and dynamical interactions. Any 
dependence of A/bh on the initial seed mass is largely erased. Dwarf galaxies undergo 
a quieter merger history, and as a result, at low masses the MBH occupation fraction 
and the distribution of MBH masses still retain some "memory" of the original seed mass 
distribution. The signatures of MBH for mation mechanisms will consequently be stronger 
in dwarf galaxies ( Volonteri et alj2008j ). Van Wassenhove et al. (2009) find that for the 
most part MBHs hosted in Milky Way satellites retain the original 'seed' mass, thus 
providing a clear indication of what the properties of the seeds were. MBHs generated 
as 'massive seeds' have larger masses, that would favour their identification, their typical 
occupation fraction is lower, being always below 40% and decreasing to less than a % 
for 'true' dwarf galaxy sizes. Population III remnants have a higher occupation fraction, 
but their masses have not grown much since formation, making their detection harder. 

This can be understood by exploring the establishment and evolution of the empirical 
correlation between black hole mass (Mbh) and velocity dispersion (a) as a function 
of cosmic time in galaxies of different mass (| Volonteri fc Nataraian 2009). The massive 
end of the A/bh — o relation is established early, and lower mass MBHs migrate onto it 
as hierarchical merging proceeds. The slope and scatter of the relation at the low-mass 
end instead retains memory of the seeding mechanism. We find that if MBH seeds are 
massive, ~ 10 5 Af Q , the low -mass end of the A/bh — cr flattens towards an asymptotic 
value, creating a characteristic 'plume'. This 'plume' consists of ungrown seeds, that 
merely continue to track the peak of the seed mass function down to late times. 



126 



Marta Volontcri 



References 

Abel T., Bryan G. L., Norman M. L., 2000, ApJ, 540, 39 

Arun K. G., Babak S., Berti E., Cornish N., Cutler C, Gair J., Hughes S. A., Iyer B. R., Lang 
R. N., Mandel I., Porter E. K., Sathyaprakash B. S., Sinha S., Sintes A. M., Trias M., Van 
Den Broeck C, Volonteri M., 2009, Classical and Quantum Gravity, 26, 094027 

Barth A. J., Ho L. C, Rutledge R. E., Sargent W. L. W., 2004, ApJ, 607, 90 

Barth A. J., Martini P., Nelson C. H., Ho L. C, 2003, ApJL, 594, L95 

Baumgarte T. W., Shapiro S. L., 1999, ApJ, 526, 941 

Begelman M. C, 2009, ArXiv e- prints 

Begelman M. C, Shlosman I., 2009, ApJ, 702, L5 

Begelman M. C, Volonteri M., Rees M. J., 2006, MNRAS, 370, 289 

Bond J. R., Arnett W. D., Carr B. J., 1984, ApJ, 280, 825 

Bromm V., Coppi P. S., Larson R. B., 1999, ApJL, 527, L5 

Bromm V., Coppi P. S., Larson R. B., 2002, ApJ, 564, 23 

Bromm V., Loeb A., 2003, ApJ, 596, 34 

Clark P. C, Glover S. C. O., Klessen R. S., 2008a, ApJ, 672, 757 
Clark P. C, Glover S. C. O., Klessen R. S., 2008b, ApJ, 672, 757 
Devecchi B., Volonteri M., 2009, ApJ, 694, 302 
Ferrarese L., Ford H., 2005, Space Science Reviews, 116, 523 
Fryer C. L., Woosley S. E., Heger A., 2001, ApJ, 550, 372 
Gaburov E., Lombardi J., Portegies Zwart S., 2009, ArXiv e-prints 

Glover S. C. O., Clark P. C, Greif T. H., Johnson J. L., Bromm V., Klessen R. S., Stacy 
A., 2008, in Hunt L. K., Madden S., Schneider R., eds, IAU Symposium Vol. 255 of IAU 
Symposium, Open questions in the study of population III star formation, pp 3-17 

Giiltekin K., Richstone D. O., Gebhardt K., Lauer T. R., Tremaine S., Aller M. C, Bender R., 
Dressier A., Faber S. M., Filippenko A. V., Green R., Ho L. C, Kormendy J., Magorrian 
J., Pinkney J., Siopis C, 2009, ApJ, 698, 198 

Hoyle F., Fowler W. A., 1963, MNRAS, 125, 169 

Kudritzki R.-P., Puis J., 2000, ARA&A, 38, 613 

Lodato G., Natarajan P., 2006, MNRAS, 371, 1813 

Madau P., Rees M. J., 2001, ApJL, 551, L27 

McKee C. F., Tan J. C, 2008, ApJ, 681, 771 

Omukai K., Schneider R., Haiman Z., 2008, ApJ, 686, 801 

Palla F., Zinnecker H., Maeder A., Meynet G., eds, 2002, Physics of star formation in galaxies 
Peterson B. M., Bentz M. C, Desroches L., Filippenko A. V., Ho L. C, Kaspi S., Laor A., Maoz 

D., Moran E. C, Pogge R. W., Quillen A. C, 2005, ApJ, 632, 799 
Rees M. J., 1978, in Berkhuijsen E. M., Wielebinski R., eds, Structure and Properties of Nearby 

Galaxies Vol. 77 of IAU Symposium, Emission from the nuclei of nearby galaxies - Evidence 

for massive black holes, pp 237-242 
Saijo M., Baumgarte T. W., Shapiro S. L., Shibata M., 2002, ApJ, 569, 349 
Santoro F., Shull J. M., 2006, ApJ, 643, 26 

Schneider R., Omukai K., Inoue A. K., Ferrara A., 2006, MNRAS, 369, 1437 
Sesana A., Volonteri M., Haardt F., 2007, MNRAS, 377, 1711 

Shapiro S. L., 2004, in Ho L. C, ed., Coevolution of Black Holes and Galaxies Formation of 

Supermassive Black Holes: Simulations in General Relativity, pp 103 — h 
Shibata M., Shapiro S. L., 2002, ApJ, 572, L39 
Shlosman I., Frank J., Begelman M. C, 1989, Nature, 338, 45 
Volonteri M., Gnedin N. Y., 2009, ApJ, 703, 2113 
Volonteri M., Haardt F., Madau P., 2003, ApJ, 582, 559 
Volonteri M., Lodato G., Natarajan P., 2008, MNRAS, 383, 1079 
Volonteri M., Natarajan P., 2009, MNRAS, pp 1476-+ 
Yoshida N., Omukai K., Hernquist L., Abel T., 2006, ApJ, 652, 6 
Yu Q., Tremaine S., 2002, MNRAS, 335, 965 



